Pattern selection of cracks in directionally drying fracture 
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Abstract 

We study pattern selection of cracks in directionally drying fractures by an- 
alyzing the experimental systems recently devised by C. Allain and L. Limat. 
[Phys. Rev. Lett. 74, 2981 (1995).] Proposing a simple picture of crack 
formation, we clarify the mechanism of how cracks array regularly and find 
that the interval between neighboring cracks is proportional to the 2/3 power 
of the cell thickness. This result explains well the experimental data of Allain 
and Limat. 
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The study of fracture has developed greatly JT| since Griffith wrote a breakthrough 
paper in 1920 0. In particular, recent progress is due largely to the development of well- 
controlled experimental systems @,f§]. Recently, Allain and Limat have studied periodically 
aligned crack patterns by devising an experimental system consisting of directionally drying 
fractures ||. In their experiment, a colloidal suspension was put into a rectangular cell 
in which one surface was left open in order to allow evaporation, and after a short time, 
periodically aligned cracks were observed ||. Similar phenomena have been observed in 
experiments on drying fracture with different geometries ||. Since a crack cannot adjust its 
position after it appears, these experimental results cannot be explained in the same way 
as periodic pattern formations in convective systems and reaction diffusion systems . We 
are thus led to consider the mechanism of crack formations in drying fracture. 

The question we address is to determine the interval between neighboring cracks, A. 
Although Allain and Limat gave a theoretical estimate of A and confirmed that its value is 
of the same order as the experimental result ||, their theory leads to a relation between A 
and the system thickness H as X/H ~ A — B In H, which does not fit their experimental data 
well. In this Letter, we propose a simple picture of crack formation in directionally drying 
fracture. Based on this picture, we clarify the reason why cracks are formed periodically 
and derive the scaling relation A ~ H 2 ^ 3 . Further, by comparison with the experimental 
data, we confirm the validity of the scaling relation. 

The experimental configuration we analyze is illustrated in Fig.l. We assume that the cell 
under consideration extends semi-infinitely in the positive y direction and has boundaries at 
y = (front surface), x = ±L/2, and z = ±H/2. Further, we are interested in the limiting 
case H <^ L. 

We first consider the water distribution in the material without cracks. Let <f> be the 
water volume fraction. We assume that water evaporates from the front surface at a rate 
{<P ~ <Poo)J, where 0oo is the equilibrium value of the water volume fraction, and that there 
is no flux at the other boundaries. In a bulk region, a diffusion current proportional to the 
gradient of <fi is assumed to arise. Then, the time evolution of <p is given by the diffusion 
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equation 



where D is a diffusion constant. As an initial condition, we assume that takes a constant 
value 0o > 0oo- Then, is independent of (x,z), and <f)(y,t) is expressed by using a Green 
function of the diffusion equation. Without knowing the explicit form of 0(y, t), we do know 
the following general aspects of its behavior which are important for the argument given 
below. First, when t is fixed and y is increased, (j>(y,t) monotonically approaches 0o with 
the characteristic variation scale \f~Dt. Second, 0(0, t) approaches 0oo with the time scale 
t e = D/J 2 . Then, the length scale reached by the diffusion due to the evaporation at the 



front surface, which is denoted by £, is estimated as £ = \jDt e = D/J. In the argument 
below, £ 3> L will be assumed so that we can concentrate on the idealized case that the 
evaporation at crack surfaces does not cause inhomogeneity in the x direction. 

We next discuss elastic properties of the material. From a macroscopic viewpoint, the 
material can be regarded as a homogeneous elastic medium. We thus apply the linear elastic 
theory to the calculation of the macroscopic stresses. Here, the stresses a xx , a yy and a zz 
are proportional to the sum of the corresponding strains and the volume shrinkage rate, C, 
due to the evaporation under the stress free boundary condition, where C is assumed to be 
given by 

^ 0-00 

C = -a— - — . (2) 

00 

We note that elastic constants and a may depend on the volume fraction of suspensions, 
which is fixed at t — 0, but is not strongly dependent on 0. Since elastic fields vary much 
faster than the diffusion field C, values of stress fields are determined adiabatically by 
the profile 0. Then, in principle, we can calculate stresses for a given 0(y, t) under given 
boundary conditions. 

Now, we consider a "macroscopic" crack formation. In the system under consideration, 
the seed of a crack is not supplied externally, but rather there are microscopic inhomo- 
geneities at the length scale of the particle radius. One may regard the inhomogeneity as 



micro-cracks embedded randomly in the elastic medium. If we assume such a picture, the 
Griffith criteria can be applied to the understanding of macroscopic crack formation [0. 
That is, a micro crack can grow if the energy release rate (per unit length) for the micro 
crack growth exceeds the surface energy (per unit length), and when the energy release 
rate increases further as the micro-crack grows, the crack grows acceleratedly and finally 
becomes macroscopically observable. This process corresponds to a macroscopic crack for- 
mation. However, due to the randomness of the positions of micro cracks, it is difficult to 
calculate precisely the energy release rate of each crack. Thus, a coarse-grained picture is 
needed. Since we believe that macroscopic crack formation does not depend on microscopic 
details and can be described by physical quantities defined on the macroscopic scale, we 
propose a following hypothesis: The relevant quantity for macroscopic crack formation is 
the sum of the elastic energy of macroscopic stress fields and the surface energy, and if the 
total energy decreases with crack formation, the crack appears macroscopically. Inspection 
of the hypothesis from a microscopic viewpoint may be an important study, but here we 
proceed to a discussion of crack formation processes based on this hypothesis. 

We discuss the time evolution of the system. First, the material tends to shrink due 
to the evaporation. However, since the displacement at the horizontal boundaries is zero, 
internal stress proportional to C are created. In the case that the material is sufficiently 
thin, we expect that a crack is formed along the ^-direction at the front surface where the 
stress takes a maximum value. The crack breaks the front surface and extends along the y 
direction until the energy release rate for the crack extension is equal to the surface energy. 
Therefore, the crack spacing observed in the experiment || is determined during the crack 
formation process at the front surface. Hereafter, we will discuss the crack formation at the 
strip corresponding to the front surface. Further, for simplicity, we assume that separations 
at the vertical boundaries (x = ±L/2) occur before crack formations in the bulk. Note 
that the system after such separation is equivalent to one under stress free conditions at the 
vertical boundaries. Whether this assumption is realistic or not depends on the conditions 
at the interfaces between the elastic material and the cell. However, as we will see later, 
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even in the case that the separation from the boundaries never occurs, the result of pattern 
selection of cracks is unchanged. 

As the evaporation proceeds, crack formation occurs first at a time t\ satisfying the 
equation 

max[£(L, 0(h)) - E(l, 0(h)) - E(L - I, 0(h))} = TH, (3) 

where T is the surface energy per unit length, and E(L,C(t)) is the elastic energy of the 
material with the horizontal length L and the volume shrinkage rate C(t) at time t. The 
value of I maximizing the quantity E(L, C(ii)) — E(l, C(ti)) — E{L — I, C(ti)), denoted by 
Li, specifies the position where the crack is formed. From the geometrical symmetry of the 
problem, we expect Li — L/2. This implies that a crack is formed first at the center of the 
strip. (This will be confirmed later.) 

In this way, at a time ti, there are two stripes with a horizontal length L/2. We should 
notice here that each strip has the same boundary conditions as the original one. Thus, by 
replacing variables (ti, L) in Eq.(|3|) with (i 2 , Li), we know the time t 2 at which the next crack 
formations occur at the center of each strip with the length Li/2. (Note that these crack 
formations occur simultaneously.) As bisected strips are further bisected in succession, this 
process repeats, producing strips with equal horizontal length until the evaporation finishes. 
This is the reason the cracks form periodically in space. The interval between cracks A, 
which is identical to the horizontal size of strips at t = oo, is determined by the maximum 
length 2~ n L (n : integer) shorter than the length A* satisfying 

E(2K,C 00 )-2E(K,C 00 ) = TH. (4) 

Also, as easily checked, A satisfies the inequality 

y < A < A,. (5) 

Therefore, we can estimate the value of A if we can succeed in deriving an expression for 
the elastic energy of the strip. In the following paragraphs, we will derive an expression for 
E(l,C). 
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In order to evaluate the expression of the elastic energy of the strip with a horizontal size 
I, we propose to consider a quasi-one dimensional spring network which is composed of a 
chain of iV springs along the center line and vertical springs connecting each node to a fixed 
position at the boundary (see Fig.0). The ends of the horizontal springs have no constraints 
because free boundary conditions are imposed at the vertical boundaries of the strip. Here, 
the natural lengths of horizontal and vertical springs at t = are given by a = l/N and 
H/2, respectively, and the attachment points are assumed to be positioned regularly, with 
a period a. Further, the spring constants of horizontal and vertical springs are denoted by 
k\ and respectively. This effective spring network resembles Meakin's model || when 
the model is supplemented with a breaking rule. We note, however, that according to our 
picture of crack formation discussed above, the horizontal spring is broken deterministically, 
as in Hayakawa's model ||, not probabilistically, as in Meakin's model. 

The elastic energy of the spring network is expressed by 
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Here, Ui is the displacement of the i-th node from the reference point x\ = ia — 1/2. Note 
that the linear shrinkage rate of springs is given by C/3. Under the assumption (ui/H) 2 
C/3 -C 1, expanding Eq.(^|) in Ui/H and ignoring terms of higher order than (ui/H) 2 , Eq.(||) 
reduces to 
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Further, in order to make it possible to develop an analytical argument, we take the contin- 
uum limit (a — > with fixing L) of Eq.(|7]). First, we introduce a variable u well-defined in 
this limit by 



Ut 



Hu 



(8) 



Then, when this expression is substituted into Eq.(^), E should not depend on a in this 
limit. By noting that Ui + \ —Ui = avl + 0(a 2 ) where the prime refers to differentiation with 



respect to the argument Xi/H, this requirement leads to the conditions k± ~ 1/a and &2 ~ a- 
Further, since it seems natural to assume that k\ and k-i do not depend on the horizontal 
length I, we can express k\ and ki in this limit as 
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where K is related to the Young modulus of the two dimensional elastic material, and k is 
a non-dimensional quantity of order unity. As a result, we obtain the following expression 
of the elastic energy in the continuum limit: 
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Here, x refers to Xi/H in the continuum limit. The equilibrium value of the displacement 



field u is obtained by minimizing E. Thus, u is given by a solution of the equation 



u" 



kC ~ n 
— u = 0, 
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(11) 



under stress free boundary conditions at both side ends (x = ±l/2H) : 

0. 



u' 
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(12) 



Solving this differential equation, we get an analytic expression for the displacement field u 



as 



C sinh(i 



cosh( v 3 2H/ 

where —1/2H < x < 1/2H. Here, note that Eq.flTS]) is valid only in the case 
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because we have assumed u 2 <^ C/3 <^ 1 in the derivation of Eq.(^). By substituting 
Eq. (p^3|) into Eq.(|l(]), an expression for the elastic energy E(l,C) is derived as 
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Further, under the condition Eg. (fR|) , Eq.([L^) becomes 
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As expected, it is easily confirmed that E(L, C) — E(l, C) — E(L — I, C) is maximum when 
I = L/2, because we find, using Eq. (p~5f) , that dE(l,C)/dl is a monotonically increasing 
function of I. Thus from Eqs.(f|)and (|16D, we obtain 

*l = lH\ (17) 

with 

losr 

7 = 7^ < 18 > 



Here, this relation is valid only when J reC/3A*/ H <C 1. Using Eq.(|5|), the length A between 
cracks periodically aligned is determined to be 

A ~ ^H 2 I\ (19) 



where Eq.(0) holds when H > H c = 12V3kC' 3/2 T / K . In Fig.|., the experimental data 
of Allain and Limat is plotted together with the theoretical curve given by Eq. (|Hj), 
where the value 7 = 5 x 10~ 4 [m] was used to best fit the experimental data. Using the 
experimental data (0o — 0)/0o ~ 0.1 in || and assuming a in Eq.@ is of order unity, we 
obtain ~ 0.1. Thus, from the fitting value of 7, Y/K is estimated as 10~ 8 [m], which 
is of the same order as the length scale of the microscopic inhomogeneity, i.e., the radius of 
colloidal particle. We believe that this result demonstrates the consistency of our theory. 
We here comment on the pattern selection of multiple cracks in heated glass plates 



recently studied by Yuse et al. |10|. In their experimental set-up, many seeds of cracks were 
prepared in the bottom of a heated glass plate which is dipped into water. Then, the spacing 
A of cracks which finally survive was measured for the thermal diffusion length d. Motivated 
by this experiment, Hayakawa carried out numerical simulations of a spring network model 
and found the scaling relations A ~ d 2 ^ 3 for the limit of large d fl 1[| . This result can easily 



be interpreted by replacing H in the model for directionally drying fracture with d. 
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Let us turn again to directionally drying fracture. In the analysis developed above, we 
postponed two problems. First, the diffusion length £ was assumed to be much larger than 
the horizontal width L. When £ is smaller than L, the effect of the evaporation at crack 
surfaces must be taken into account. In such a case, the above simple picture cannot be 
applied to the crack formation process because inhomogeneity of the water volume fraction 
arises. In this case, the crack spacing may have some distribution. Still, we believe that 
the scaling relation given by Eq.(|l9]) will give a first order approximation for the averaged 
value of crack spacings. Second, we assumed that the separation of the strip from the 
vertical boundary occurs before the appearance of cracks. However, even in the case that 
the separation from the boundaries never occurs, the vertical boundaries of the slice of 
material lying between cracks after their formation will satisfy the stress free boundary 
conditions. Therefore, by applying the above discussion to such a slice, we can again obtain 
the result Eq.(0). 

Finally, we briefly discuss the crack patterns which appear subsequent to the crack 
formation. As soon as a crack is formed, it extends along the y direction and is arrested 
at the position where the energy release rate for the crack is equal to the surface energy. 
After an evaporation time D/J 2 , there is no crack formation at the front surface. Then, one 
may expect that the individual crack tips will all extend to the same position. However, 



as shown by H-A. Bahr et al. [0, a zig-zag type or more complicated arrangement occurs 
when a certain condition is satisfied. In fact, we can see such patterns in the experiment 
performed by Allain and Limat ||. Further, in their experiment, a more interesting crack 
pattern appeared due to secondary branching from the elongating cracks. Such patterns will 
be studied by devising an extended version of the present model. 
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FIGURES 

FIG. 1. Schematic view of experimental set-up. The axes x, y, z are defined as seen. The drying 
process takes place only at the front surface, whose dimension is defined as L x H. Cracks originate 
at the surface and extend in the y direction. 

FIG. 2. Schematic view of effective one dimensional model for drying fracture describing frac- 
ture pattern formation at the front surface (see Fig|l|). Each component (circle) is connected with 
its nearest neighbor components by k\ springs and with upper and lower plates by hi springs. In 
the drying process, k\ springs may break. 

FIG. 3. Interval of periodic patterns in drying fracture (vertical axis) vs. height of strips 
(horizontal axis). The solid line ( r y 1 / 3 H 2 ^ 3 ) is the theoretical curve, where 7 1 / 3 = 8.2 is determined 
to fit the experimental data. Filled circles with error bars are experimental data taken from Allain 
and Limat ||. The unit length for each axes is fim. 
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T.S.Komatsu and S.Sasa Figure 1. 




T.S.Komatsu and S.Sasa Figure 2. 



